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Abstract: We study the response of the Sakai-Sugimoto holographic model of large N c QCD 
at nonzero temperature to external electric and magnetic fields. In the electric case we find a 
first-order insulator-conductor transition in both the confining and deconfining phases of the 
model. In the deconfining phase the conductor is described by the parallel 8-brane-anti-8- 
brane embedding with a current of quarks and anti-quarks. We compute the conductivity and 
show that it agrees precisely with a computation using the Kubo formula. In the confining 
phase we propose a new kind of 8-brane embedding, corresponding to a baryonic conductor. In 
the magnetic field case we show that the critical temperature for chiral-symmetry restoration 
in the deconfined phase increases with the field and approaches a finite value in the limit of an 
infinite magnetic field. We also illustrate the nonlinear behavior of the electric and magnetic 
susceptibilities in the different phases. 
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1. Introduction and Summary 

Holographic descriptions of large N c QCD-like theories have received a great deal of attention 
recently (for a recent review see [1]). The models involve the physics of Nf flavor-brane 
probes in near-horizon geometries of N c color-branes. The Sakai-Sugimoto D4-D8-D8 model 
[2] in particular is very attractive in that it has a simple geometric description of chiral 
symmetry breaking, and appears to be rich enough to incorporate, at least qualitatively, all 
the low energy features of QCD. The properties of the mesons [2, 3] and baryons [4], the 
resolution of the U{\)a problem [5], and the phase diagram at nonzero temperature [6, 7], 
nonzero baryon chemical potential [8, 9], and nonzero isospin chemical potential [10], all 
exhibit many similarities with QCD. While the original model had only massless quarks, it 
can be generalized to nonzero quark mass [11]. 

The Sakai-Sugimoto model does not include a true electromagnetic gauge field, but we 
can mimic the effect of one using the Abelian part of the flavor symmetry. In this paper we will 
analyze the Sakai-Sugimoto model at finite temperature with an Abelian background gauge 
field in the diagonal U(l)v part of the D8-D8 gauge group. This gauge field is holographically 
dual to the baryon number current of the four-dimensional gauge theory. Except for the one- 
flavor case, this is not the same as an electromagnetic field, since "up" and "down" type 
quarks have the same charge. Nevertheless, we expect the physics to be qualitatively similar 
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to that of a background electromagnetic field. Background flavor gauge fields have been 
studied previously in the ftf = 2 supersymmetric theory corresponding to 7-brane probes in 
the 3-brane background [12, 13, 14]. 

We would like to determine the effect of the electric and magnetic fields on the phase 
diagram of the theory and compute the response coefficients (susceptibilities or conductivities) 
of the different phases. We will do this by analyzing the 8-brane embeddings in the presence 
of the appropriate background world-volume gauge field. The value of the 8-brane action is 
identified in the usual way with the appropriate thermodynamic potential. 

In the case of an electric field we find an insulator-conductor type transition, which 
in the deconfined phase generalizes the chrial symmetry breaking-restoration transition at 
zero field. The critical temperature decreases with the field. Since we use the DBI action, 
there is a temperature-dependent maximal value of the electric field. However, the transition 
to the conductor occurs at a smaller value of the field. We compute the conductivity in 
the conducting phase and show that it agrees precisely, in the zero field limit, with what is 
expected from the Kubo formula. This agreement extends to the nonzero density case as well. 
Somewhat surprisingly, we also find an insulator-conductor transition in the confined phase, 
where there is no chiral-symmetric phase. We propose that the conducting phase corresponds 
to an embedding in which the 8-brane and anti-8-brane are geodesically parallel and connect 
at a cusp at the "tip of the cigar". We compute the conductivity of this embedding and 
determine the phase diagram. The current in this case is carried by baryons and anti-baryons. 

With a magnetic field we observe that the critical temperature for chiral symmetry 
restoration increases with the field, in agreement with expectations for QCD. We find that 
the critical temperature approaches a finite value in the limit of infinite field, which differs 
from the behavior in the J\f = 2 theory found in [13, 14]. 

The paper is organized as follows: In section 2 we review briefly the Sakai-Sugimoto 
model and set up our conventions; sections 3 and 4 deal with, respectively, the physics of 
background electric and magnetic fields. 



2. Review of the Sakai-Sugimoto model 

The model consists of N c D4-branes in Type IIA string theory wrapping a circle with anti- 
periodic boundary conditions for fermions, Nf D8-branes at a point on the circle, and Nf 
anti-D8-branes at another point on the circle. At energies well below the Kaluza-Klein scale 
the spectrum on the D4-branes is precisely that of massless four-dimensional "QCD" , with 
N c colors of gluons and Nf flavors of quarks. The holographic dual description comes about 
by taking the large N c limit, in which the D4-branes are replaced by their near-horizon 
supergravity background. It is most convenient to work in units in which the curvature radius 
is fixed to unity, R = (-Kg s N c ) 1 ^ \fa' = 1. Taking x\ as the circle direction, X4 ~ £4 + 27ri?4, 
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the background is given by 



3 

U 1 



g s u 



dx 2 + dx 2 + f(u)dxl) + 
!( F 4 = 37r(a / ) 3/2 /V c dQ 4 , 



_ 3 

u 2 



du 2 



+ u 2 dnj 



(2.1) 



where /(it) = 1 — (u KK /u 3 ), and Uft-^; = 4/(9i? 4 ). The (14,2:4) subspace is topologically a 
cigar (or disk), with a tip at u = ukk- In the dual gauge theory this ultimately implies that 
the gluons are confined. The gravitational description is valid as long as A>i?4, where A is 
the five-dimensional 't Hooft coupling, which in our units is given by 

4 



A = 4Tig s N c Va' 



a' 



(2.2) 



The 8-branes and anti-8-branes are treated as probes in this background, and their em- 
bedding determines the flavor physics in the dual gauge theory. Due to the topology of the 
background the 8-branes and anti-8-branes must connect into a smooth U-shaped configura- 
tion at some radial position uq > Ukk (fig- l a )- This reflects the spontaneous breaking of 
the U(Nf)n x U{Nf)L chiral symmetry to the diagonal U(NAy. The form of the embedding 
can be determined from the DBI action of the 8-branes in this background, 



d 4 x 



du u 



S = M 

The normalization constant is given by 

TV = 2N f V^n 4 T 8 



f(u)(x' 4 (u)) 2 + 



1 



u 3 f(u) 



1/2 



4N f Vi 



(2.3) 



(2.4) 



3(2vr)6(a') 9 / 2 9s 

where V4 is the volume of 4d spacetime, Q4 is the volume of a unit 4-sphere, and Ts is the 
8-brane tension. The factor of 2 corresponds to the two halves of the embedding (8-branes 
and anti-8-branes) along u. The equation of motion for the embedding 2:4(14) gives 

1-1/2 

x' 4 (u) 



1 



Therefore at large u 



0:4(14) 



u s f(u) 



Uof(u ) 



1 



L 
2 



9 u 9 / 2 ' 

where L is the asymptotic 8-brane-anti-8-brane separation, 



L 



du 0:4(14) , 



(2.5) 



(2.6) 



(2.7) 



and c is the constant of the motion associated with 2:4(14). This constant corresponds to the 
curvature of the 8-brane, and is related to uq by 



it, 



uq) 



(2.8) 
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At nonzero temperature there are two possible backgrounds. For T < l/(2-7ri?4) the 
dominant background is the Euclidean continuation of (2.1) with Xq ~ Xq + 1/T. In this 
background the 8-brane embedding is the same as above. The glue sector is therefore confined, 
and chiral symmetry is broken. For T > l/(27ri?4) the dominant background is given by (2.1), 
with the roles of x^ and Xq exchanged and with f(u) = 1 — (uj,/u 3 ), where ut = (47rT/3) 2 . 
Here the (u, X4) subspace is topologically a cylinder, with a horizon at u = ut, which in the 
dual gauge theory implies deconfinement. In this background there are two possible 8-brane 
embeddings: a U-shaped embedding (fig. lb), that satisfies 



x' 4 (u) 



u 



3/2^ 



i4f{u Q ) 



(2.9) 



and a parallel 8-brane-anti-8-brane embedding with x' 4 (u) = (fig. lc). For T < 0.154/L the 
U embedding dominates and therefore chiral symmetry is broken, but when T > 0.154/L the 
parallel embedding dominates and chiral symmetry is restored. The intermediate phase of 
deconfinement and chiral symmetry breaking appears only when this critical temperature is 
higher than the deconfinement temperature l/(2-7ri?4), namely when L < 0.97 R4. Both the 
confinement /deconfinement and chiral symmetry breaking/restoration transitions are first 
order. 




a b c 

Figure 1: 8-brane embeddings and phases of the Sakai-Sugimoto model: (a) confined, 
broken chiral symmetry (b) deconfined, broken chiral symmetry (c) deconfined, restored 
chiral symmetry 



The 8-brane curvature c is an order parameter for the chiral symmetry transition in the 
deconfined phase: it vanishes in the chiral-symmetric parallel embedding, and is given by 
(2.8) in the chiral symmetry breaking U embedding. 1 By studying the dependence of L on c 
one finds a turn-around behavior typical of a first order phase transition. This is easily seen 
from the asymptotic behavior of L at small and large c (see fig. 2 for a numerical plot of L 

lr This is not the usual chiral symmetry order parameter. The 8-brane curvature c is related to the expec- 
tation value of a quark quadra-linear operator [15]. The usual chiral symmetry order parameter is the quark 
bi-linear (qq), which is given by the normalizable mode of the flavor bi-fundamental scalar field [11]. 
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vs. c at a fixed temperature, and also [7]), 



c ^0 (u ^u T ) : L(c,T)~c/T 9 
c — ► oo (uq — > oo) : L(c, T) ~ c -1 / 8 . 

This implies that there is a maximal value of the asymptotic separation L max , which depends 
on T. For L < L max there are actually two U embedding solutions, and for L > L max there 
are none. Alternatively, since L is a monotonically decreasing function of the temperature, 
there is a maximal temperature T max at any fixed L. For T < T max there are three solutions in 
all: the parallel embedding, and the two U embeddings. Evidently, one of the U embeddings 
must be an unstable solution. At T = T max the unstable U embedding merges with the stable 
one, and at higher temperatures the two disappear, leaving only the parallel embedding. The 
transition from the stable U embedding to the parallel embedding occurs at a temperature 
lower than T max . 

L 

1 

.75 
. 5 
0.25 

0.025 0.05 0.075 0.1 C 

Figure 2: L vs. c for the U embedding in the deconfined phase (T = 0.14) 




3. Electric Field 

In this section we will study the response of the model to an external electric field E, by 
turning on an appropriate background value for the Abelian gauge field component of the 
unbroken U{Nf)y gauge field in the 8-brane worldvolume. We normalize this field as follows 

A=±~TrA, (3.1) 

where A is the U(Nf)y gauge field. 2 Anticipating a current in the direction of the background 
field, we make the ansatz (in Euclidean space) 

Ai(x ,u) = -iExfi + F(u), (3.2) 

where the u dependence encodes the current in the usual holographic fashion [12]. 

2 The canonical normalization is A — y/2/Nf TrA, but then the quarks carry 1/ y/2Nf units of charge. 
In our normalization the kinetic term has an extra factor of 2Nf, but the quarks carry unit charge. This 
normalization is more suitable to mimic an electromagnetic field. 
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3.1 Deconfined phase 



We begin in the deconfining background, which dominates when T > l/ftirR^). The DBI 
action for the 8-branes with the gauge field in (3.2) is given by 



S = M I duu 



f(u)(x 



+ 



1 



f{u)u c 



+ 



l\2 



(3.3) 



where we have defined the dimensionless quantities a% = 2ira' ' A\ and e = 27ra'E. The 
asymptotic behavior of the gauge field is given by 



ai(xo, u) ~ —iex\ 



2 j 



3 u 3 / 2 ' 



(3.4) 



where j is the conserved charge associated with a\, namely the baryon number current. The 
physical dimensionful current is given in our units by J = (2irot!j\f/V<£)j. In terms of the 
current the action becomes 



S = N I duu A J ( f{u)x' A 2 + 



1 



ur 



/(«) 



ur 



/(«) 



r 



a'- 



(3.5) 



This form of the action displays a generalization of the usual limiting electric field of the 
flat space DBI action. For a vanishing current, the action is complex, and the embedding is 
unphysical, if the second factor becomes negative anywhere in the integration region. However 
this can be fixed by turning on a current such that the third factor changes sign at the same 
point. This gives an equation relating the current and the electric field [12]. 

As in the zero field case there are two types of embedding. Consider first a U embedding 
with a vanishing current, j = 0. The solution satisfies 



x' A (u) 



1 



u 8 /(«) 



e 



«5 /(«o) 



e- 



-1/2 



(3.6) 



and the large u behavior is given by (2.6), with 



c = u o\ f( u o) 



(3.7) 



The action of this solution is given by 

duu 5 /\ll 



u 3 f(u) 



4 ( /(no) - $ o 



u 8 (f(u) 



-1/2 



(3.8) 



The U embedding with j = is always physical: since c is real, the solution satisfies e 2 < 
Uq/(uo), and the action is real. A current may be turned on, as long as j 2 < ^/(ito), but 
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this increases the action, so the dominant U embedding has j = 0. In the gauge theory this 
embedding therefore corresponds to a chiral-symmetry breaking, insulating phase. 

In fact the U embedding satisfies an even tighter bound on the electric field than above. 
As in the zero field case, at fixed values of e and T, there is a maximal value of L as a function 
of c for the U embedding. Since L is a monotonically decreasing function of e, this implies a 
maximal value of e as a function of c at fixed values of L and T. The maximal value e max is 
attained at some c > 0, which, using (3.7), implies that e^ ax < Uq/(uo)- For e > e max there 
are no U embedding solutions, and we therefore expect a phase transition to occur, at fixed 
T and L, at some value of c smaller than G max ■ 

In the parallel embedding x'^{u) = 0, and the action is given by 



roc 

SW=N duu^ 2 

J Ut 



f(u) - 4 

T • (3-9) 



The numerator is negative for u 3 < e 2 + which is always in the range of integration. The 
only way to ensure a real action in this case is for the denominator to become negative at the 
same u. This requires a nonzero current given by 

j = e ( e 2 + 4 )i/3 ( 310 ) 

The parallel embedding therefore describes a chiral-symmetric conducting phase in the gauge 
theory, and the conductivity is given by 

' = i = t (e2 + * 1/3 = *%nr- f 1 + ^ 1/3 • (3 - n) 

where we have defined a new dimensionless variable e by 

e 27 E 

e = ^[ 2 ~ stoats ' (3 - 12) 

To determine which phase dominates as a function of the temperature and electric field 
we should in principle compare the electric free energies of the two solutions, which are in 
turn defined by the Euclidean 8-brane action of the solutions 

T e {L,e,T) =TS[x i (u),a 1 (u),T] EO M ■ (3.13) 

However this is not quite right for the parallel embedding. First of all, the conducting phase 
is not in equilibrium. There is a steady state current with a finite conductivity, meaning 
that energy must be constantly added to the system. This energy is dissipated into the 
gluon "bath", which, in general, raises the temperature. At large N c , however, this effect 
is negligible. While the dissipated energy is 0(N C ), there are 0(N 2 ) gluons among which 
to distribute this energy, so the temperature rise is only 0(N~ l ). But even ignoring the 
dissipation, one still needs to subtract the kinetic energy of the current carriers, which should 
not be taken as part of the budget at the phase transition. 
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Alternatively, we can get around this problem using a Maxwell-like construction for the 
order parameter 



dL 



(3.14) 



In the parallel embedding c = for any L. In the U embedding the dependence of c on L can 
be determined numerically using (2.7) and (3.7). Note that here we need the full solution, 
not just the asymptotic behavior. The result is qualitatively the same as in the zero field 
case (fig. 2). The transition occurs at the value of L such that the two bounded areas are 
equal (see fig. 3 for an illustration). We can then construct the phase diagram in the (T, e) 



Figure 3: Illustration of the Maxwell construction 



plane with fixed L by repeating this procedure for various values of e and T, and finding 
the points which have the same critical L. The result, shown in fig. 4, shows a first-order 
insulator-conductor transition at nonzero temperature and background electric field. At zero 
electric field this reduces to the chiral-symmetry breaking-restoration transition. 



0.15 



o.l 



0.05 



Parallel embedding 




0.05 
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0.15 



Figure 4: Phase diagram at nonzero temperature and electric field in the deconfining 
background (L = 1) 
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3.1.1 Conductivity at finite density 

The conductivity computed in (3.11) is for the vacuum at a finite temperature. It has two 
contributions corresponding to quantum and thermal pair-creation of quarks and anti-quarks. 
The introduction of a finite charge density should also contribute to the conductivity. We 
can compute this by generalizing our analysis in the parallel embedding to the finite density 
case. At nonzero density there is a non-trivial time-component of the gauge field ao(u), and 
the associated constant of the motion is the dimensionless baryon number charge density d. 
The dimensionful charge density is D = {2-na'N /V^)d. The 8-brane action is 3 



S = M I duii 



1 



1 



+ 



K)2 



f(u)u 3 / 

In terms of the current and density, the action of the parallel embedding becomes 



(3.15) 



/(«) 



\ /(«) 



(3.16) 



The condition for reality now gives 



(e 2 + 4) 2/3 + 



( / 2 



e 2 + Un 



1/2 



(3.17) 



and the conductivity is therefore 

(2ro')W 

N f N c XT 2 



a 



+ 4) 2/3 + 



27vr 



fl + e 2 ) 2/3 + 



d 2 



,2 
d 2 



1/2 



e 2 + Uj, 

1/2 



1 + e 2 



(3.18) 



where we have defined a new dimensionless (with appropriate factors of R, which in our units 
is 1) variable d, 



d 



d 



729 



D 



u 5/2 SvNfNc A 2 T 5 ' 



(3.19) 



At zero charge density this reduces to the vacuum result (3.11). 
3.2 Conductivity and the Kubo formula 

The electrical conductivity is an example of a transport coefficient that describes the response 
of a thermodynamic system to a disturbance which takes it out of equilibrium. For a small 



3 There is also CS term of the form J 0,00,1(82(13 — 8302). However this can be consistently set zero by the 
equations of motion since ao and ai are assumed to be independent of X2 and X3. 
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disturbance, i.e. near equilibrium, transport coefficients can be related to real-time correlation 
functions at equlibrium via Kubo formulas. The electrical conductivity near equilibrium is 
related to the current-current correlator [16], 



a= Km o —TrC<(k = \k\), (3.20) 

where 



C<(k) = J d'xe-^iJ^Mx)) . (3.21) 

The correlator can in turn be computed at strong coupling using the Lorentzian AdS/CFT 
prescription of [17]. 4 It is therefore interesting to compare the result of this computation with 
the direct computation of the conductivity in the previous section. Since the Kubo formula 
gives the conductivity near equilibrium, one should compare the result with the zero electric 
field limit of (3.18). 

The current-current correlator for light-like momenta has been analyzed in the Sakai- 
Sugimoto model at both finite temperature and density in [18]. In particular the high- 
temperature, or equivalently low frequency, behavior was found to be 5 



lim^ C <(fa = |k|) = 2 < 2JV ^ AT3 V^ 

fco^O 277T 



(3.22) 



leading to a conductivity 



N f N c \T 2 I ^ 

<y= J 277r Vl + D 2 , (3.23) 

in perfect agreement with the zero field limit of (3.18). 

3.3 Confined phase 

For T < l/(2-7r .R4) the confining background dominates, and the 8-brane action with the 
background gauge field is 



8-mJ iuu'Mmw + J ±^(i-^ + w> 2 



u 3 



= » I d »»V ('<"W 2 + t4?) - 5) - 5) " ■ (3 - 24) 

At zero field the only possible embedding in the confined phase was the U embedding. However 
as we increase the electric field we encounter a puzzle. As in the deconfined phase, the U 



4 This prescription actually yields the retarded correlator C^„(k) — i J d^xe tk ' x 8(xo)([Jn(0), J v (x)\), which 
is related to the one above by C< v {k) = 21m C^(fc)/(e _fc ° /T - 1). 

5 Our definition of A is different than in [18]: Xhere = 47rA t h ere . Our D is exactly their C. 
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embedding solution at a fixed L ceases to exist above a certain value of the field. To see this 
we again study the behavior of L as a function of c. The U embedding with j = satisfies 



x' 4 (u) 



i 3 / 2 /(? 



u*f{u) 1 



_ M o/K) (i - % 



and 



-1/2 



(3.25) 



'/(«o) 1 



(3.26) 



There are two cases to consider. For e 2 < u^ K , L(c) decreases monotonically from L(0) = 
7ri?4 (the anti-podal embedding) to zero as c — > oo (fig. 5a). For e 2 > the asymptotic 
behavior of L(c) is the same as in the deconfined phase (fig. 5b), implying a maximal L for 
a given e, or alternatively a maximal e for a given L. This maximal value therefore satisfies 



< e 



. < Uq. The U embedding exists only for e < e max , and its action is given by 



5/2 



S u =N du^= 

Ju V /( 



\ "g/(m)(l-§ 

u 8 f(u 



-1/2 



(3.27) 
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Figure 5: L vs. c in the confined phase (a) e < (b) e > 



3/2 



The question is what happens when e > e ma z? There must exist a second embedding 
that takes over before reaching the maximal electric field. In the deconfined background this 
was the current-carrying parallel embedding, and we observed a first-order phase transition at 
e < e m ax- We therefore propose a new kind of 8-brane embedding in the confining background, 
which is analogous to the parallel embedding in the deconfining background. The 8-brane 
and anti-8-brane follow parallel radial geodesies and connect at u = ukk (fig- 6). In this 
"V-shaped" embedding x' 4 (n) = (and therefore c = 0) except at the tip, where there is a 



cusp (unless we are in the anti-podal embedding, which is smooth). Away from ukk this is 
clearly a solution. Its action is given by 



roo 5/2 

S V =M du^=. 

J U KK Vf(u) \ 



ir' 



(3.28) 



It follows from reality of the action that if e 2 > v? KK there must be a current given by j = e 
This embedding is therefore a conductor, with a conductivity 



5/3 



(2W) 2 AA 2/3 _ N f N c 1/3 2/3 

V 4 " 127TV3 



(3.29) 



Figure 6: The conducting V embedding in the confining background 



This proposal raises two questions. The first has to do with stability of the cusp sin- 
gularity at the tip, and the second with the identification of the current carriers. In the 
deconfined parallel embedding the current is carried by fundamental strings, namely quarks 
and anti-quarks. The only charged objects in the confined phase are baryons, so the current 
can only be carried by baryons and anti-baryons, namely by 4-branes (and anti-4-branes) 
wrapped on the S 4 . Indeed, this is precisely what we need to source the 5d magnetic field 
a'^u) dual to the boundary current j. This is similar to the situation at nonzero density 
[8], where a uniform distribution of 4-branes was used to source the 5d electric field a' (u) 
dual to the baryon number density. In that case the 4-branes created a cusp in the 8-brane 
embedding. In the present case the current corresponds to a uniform distribution of 4-branes 
and anti-4-branes (since the total charge density vanishes), located at the cusp, and moving 
at a constant velocity along x\. To understand whether this configuration is stable we need 
to understand the forces at the cusp. In the nonzero density case the upward force due to 
the 8-brane was balanced against the downward force of the 4-branes. In this case however 
the 4-branes are at the bottom of the space, and it is not clear how the balance comes about. 
Higher derivative corrections may also be relevant. We leave this as an open question. 

Assuming this is a valid solution, we can construct the phase diagram using the same 
method as in the deconfined phase. In this case we can either fix ukk and vary L and e, or 
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fix L and vary ukk and e. The two results are shown in figure 7. From the first diagram we 
see that for e 2 < u% K the U embedding dominates at all values of L, and that for e 2 > v? KK 
there is a first-order transition to the conducting V embedding at a critical L that starts at 
7ri?4 (the anti-podal embedding) and decreases with e. 




Figure 7: Electric phase diagram in the confined phase (a) fixed ukk = 0.5 (b) fixed L = 1. 
In (a) the dashed lines illustrate that for e = uJ? K ~ 0.35 the critical L = ttR^ ~ 3. 



The second diagram also provides some interesting insight. First, we note that, for a 
fixed L, the critical electric field in the limit ukk — > in fig. 7b is exactly the same as the 
critical field in the deconfined phase in the limit T — ► (as seen in fig. 4). This is because 
both cases are then essentially equivalent to the non-compact model at zero temperature. 
Second, by comparing the phase diagrams in the confined and deconfined phases at the 
confinement /deconfinement transition ukk = ut, it can be seen that if we start in the 
confined conducting phase (V embedding) and raise the temperature above the deconfinement 
temperature we end up in the deconfined conducting phase (parallel embedding). In this 
transition the conductivity jumps from (3.29) to (3.11). 

3.4 Electric susceptibility 

In both the confined and deconfined phases there is an insulating phase (U embedding) in 
which the current vanishes. These phases exhibit an electric polarization and susceptibility, 
which are defined thermodynamically by 

P = -^,Xe = -^- (3.30) 

In the holographic prescription the free energy is divergent, and requires a counterterm pro- 
portional to e 2 . The susceptibility therefore requires a counterterm independent of e. A 
physically motivated scheme is to require the susceptibility to vanish at zero field, in other 
words 



Xe[e) ~ 9e 2 + de 2 



(3.31) 
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This is a measure of the nonlinearity of the vacuum in this model. The results are presented 
in figure 8. 




Figure 8: Electric susceptibility in the insulating phase of the (a) confined phase, (b) 
deconfined phase (the dashed line shows the critical field for the insulator /conductor 
transition). 



4. Magnetic field 

For an external magnetic field H our ansatz is simply 

A2(xi) = Hxi . 



(4.1) 



There will be no current and therefore no u dependence. Our main objective here is to 
determine the effect of the background magnetic field on the critical temperature for chiral- 
symmetry restoration in the deconfined phase. We will also compute the magnetic suscepti- 
bility of the vacuum in both the confined and deconfined phases. 



4.1 Deconfined phase 

We start again in the deconfined phase. The 8-brane action in the magnetic field background 
is given by 



M 4 l/(/(«K 2 + ^ 



1 + 



ft 2 



If 



where h = 2na'H. The chiral symmetry breaking U embedding now satisfies 



x' 4 (u) 



« 8 (/(u) + h - 
h 2 

3 




< /(«o) + 5 



-1/2 



(4.2) 



(4.3) 
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and has an action 



JUQ V 



h 2 



u 3 f(u) 



o /(«o) + 5 



u8 (/(«) + S) 



-1/2 



(4.4) 



In the chiral-symmetric parallel embedding ^(u) = 0, and the action is 

/•oo 

S^=Af *Wu 3 + /i 2 . 



(4.5) 



The actions of the embeddings, which define the magnetic free energies J- m (L,T,h), are 
divergent, but the difference is finite. The resulting phase diagram (for a fixed value of L) is 
shown in figure 9a. Note that unlike the electric field case, here both embeddings describe 
equilibrium states, so we can compare the actions directly. 

We observe that the temperature at which chiral symmetry is restored increases with 
the background magnetic field, and approaches a finite value in the limit of an infinite field. 
This means that above some nonzero temperature chiral symmetry is always restored. A 
similar increase in the critical temperature for the phase transition in the D3-D7 model was 
observed in [13, 14], but there is a crucial difference with our result. In the D3-D7 model the 
critical temperature diverges at a finite value of the magnetic field, which means that there 
is no phase transition for magnetic fields larger than this value (fig. 9b). It is ammusing to 
speculate whether this can be tested in real QCD, either experimentally or on the lattice. 
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Figure 9: (a) Our phase diagram with a magnetic field (b) analogous phase diagram in 
the D3-D7 model (reprinted from [14] with the authors' permission). 



4.2 Confined phase 

In the confined phase the 8-brane action is 



h 2 
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The U embedding is basically the same as in the electric case, (3.25) and (3.26), with e 2 
replaced by — h 2 . However, the sign difference guarantees that this solution is the only one 
and that it exists for all values of L and h. The action of the solution can also be read off 
from the electric case (3.27) with the above replacement, 



5/2 



S u = N / du^= 



1 + 



h 2 



u>/(«) (l + S 



-1/2 



(4.7) 



4.3 Magnetic susceptibility 

The magnetization and magnetic susceptibility are defined thermodynamically by 

BT B 2 T 

UJ 777, U J 777, 



m 



dh 



X;, 



dh 2 ' 



(4.8) 



We can compute the deviation from linearity by regularizing the susceptibility as in the 
electric case, 



Xm(h) 



+ 



dh 2 dh 2 

The results for all three phases are shown in figure 10. 



h=0 



(4.9) 
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Figure 10: Magnetic susceptibility: (a) confined, broken chiral symmetry (b) deconfined, 
broken chiral symmetry (c) deconfined, restored chiral symmetry 
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